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Abstract 

It is known that silicon is an indirect band gap material, reducing its efficiency in photovoltaic 
applications. Using surface plasmons in metallic nanoparticles embedded in a solar cell has recently 
been proposed as a way to increase the efficiency of thin film silicon solar cells. The dipole mode 
that dominates the plasmons in small particles produces an electric field having Fourier components 
with all wave numbers. In this work, we show that such a field creates electron-hole-pairs without 
phonon assistance, and discuss the importance of this effect compared to radiation from the particle 
and losses due to heating. 



I. INTRODUCTION 



Present day solar cell industry is completely dominated by the use of silicon as the 
active material. It's main advantages are availability, disposability and several decades of 
industrial metallurgical development, compared to the poisonous or rare elements of, e.g., 
GaAs. However, silicon is not an ideal material for solar cells. One disadvantage is that it 
has an indirect band gap. This means that photons with energy close to the band gap can 
only be absorbed in phonon-assisted processes. Therefore, the absorption of these photons is 
weak, and the silicon wafer can not be made too thin if one is to absorb this part of the solar 
spectrum. The material costs and limited production capacity for solar grade silicon mean 
that the thickness required in todays first generation solar cells is a significant obstacle to 
their commercial success. Also, due to the limited lifetime of the electron-hole pairs, thicker 
cells may suffer from larger recombination rate and reduced efficiency. 

The question is thus how to increase the optical path lengths of near band gap photons 
inside the silicon, without increasing wafer thickness. Several approaches have been tried, 
including texturing of the wafer front or rear surface in various patterns and on different 
length scales. For length scales larger than the wave length, the incoming light is refracted 
into angles more parallel to the wafer. ^ For length scales close to or smaller than the wave- 
length, diffraction may couple light into guided modes. ^ However, the texturing will often 
lead to an increase in surface defect states and thereby increase recombination rates. 

As an alternative to texturing, it has been proposed to place metallic nanoparticles near 
the surface of the wafer.- The nanoparticles scatter the incoming light through a surface 
plasmon resonance. Surface plasmons, or surface plasmon polaritons, are electron density 
fluctuations at the interface between a metal and a dielectric material. For a good intro- 
duction, see, e.g., Raether^ for plasmons in general, and Bohren&Huffman^ for plasmons on 
small particles. On the surface of nanoparticles, the plasmons can be excited by an incoming 
plane wave, and they exhibit a marked, tunable resonance. For frequencies near the reso- 
nance, nanoparticles have an optical cross section much larger than their geometrical cross 
section. If this resonance could be tuned to match the band gap of silicon, near bandgap 
photons could be absorbed into the plasmon state with high probability, while higher en- 
ergy photons would be unaffected. Certain progress has already been made in the plasmon 
tuning, e.g., at the University of New South Wales,-- but further development is required. 
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We believe that the energy of the surface plasmons can then be used to create electron- 
hole pairs in two ways. First, the energy can be emitted as light in directions along the 
wafer. This gives a longer optical path inside the wafer, and thereby increases the indirect 
absorption.- 

Second, the near field of the nanoparticles can excite electron-hole pairs without phonon 
assistance, the momentum being transferred to the nanoparticle. This second process has to 
our knowledge not been considered in the literature, and is the subject of the present paper. 
Our results indicate that this mechanism will give an extra contribution to the electron-hole 
pair generation, compared to estimates that only take into account the re-radiation of power, 
increasing the relative benefit of introducing the nanoparticles. 

Some of the plasmon energy goes into heating of the nanoparticles and is obviously lost. 
This loss should be compared with the losses due to the limited optical path when not 
exploiting the plasmons, or, if texturing is used to increase the optical path length, with 
increased losses due to recombination at interfaces. The plasmons will give an improved 
efficiency if the losses to heating are smaller than previous losses due to optical path length 
or recombinations. If the resonance is properly tuned, the only photons significantly affected 
will be those that would otherwise be lost. Any fraction of this near bandgap light that can 
be used efficiently contributes to a net gain for the cell. 

For particles larger than the wavelength, a large fraction of the light will be reflected rather 
than excite plasmons.- We therefore consider only particles smaller than the wavelength of 
the incoming light. For such particles, the plasmons can be approximated by a dipole mode, 
corresponding to uniform polarization of the nanoparticle.- While the dipole approximation 
is usually only accepted for particles with diameter less than one tenth of a wavelength, we 
accept it as a flrst approximation for our order of magnitude estimates. We are not aware 
of any studies of how the near fleld is changed by an interface between the layer embeddig 
nanoparticles and the active layer of the solar cell. For simplicity, we will therefore restrict 
the further discussion to the case of an electric dipole located inside an inflnite medium 
consisting of silicon. The dipole is excited by an incoming plane wave. 

The far-fleld energy radiated from the dipole represents the maximum energy that can 
be absorbed by indirect absorption. In real applications, some of this light will inevitably 
be lost. The presence of an interface may also increase the total emission,-i£ but for the sake 
of our order of magnitude estimates, we will ignore this effect. 
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The goal of this paper is to demonstrate that the direct absorption effect should be con- 
sidered when modeling the effect of plasmons, and that it may have important implications 
for the optimal sizing and positioning of the plasmons. The fact that plasmons can lead 
to an increase in efficiency has been experimentally verified.-'^ We therefore focus on the 
relative importance of the two mechanisms that could contribute to the increase, and how 
this could influence cell design considerations. 



II. THEORY 



We use classical electrodynamics to describe both the nanoparticles and the fields. The 
interaction with the silicon is described by perturbation theory, and we use the tight binding 
model and the parabolic approximation of the band gap extrema for the wave function of 
the silicon. 

We consider the incoming light to have a frequency close to the band gap of silicon. 
This corresponds to hu = 1.1 eV, or u ^ 10^^ s~^. We then get for the wavelength of 
this radiation A = 2ttc/u! ~ 1 /im , and from c = u/kp we get the photon wavenumber 
kp^ 6 ■ 10^ m~^. 

The vector potential due to a dipole is given as:— 

A = ^e-(^^'^--*)po = .AoA , (1) 

_ UJfiokpPo _ 

^° " 47r ~ kpr"" 

where po is the dipole moment, and Po = |Po|- '"is the distance from the dipole. Using 
the previous rough estimates for u and kp we get that Aq/pq ^ 10^^ Js/C^m^, while A,, is a 
dimensionless function containing all spatial dependencies of A. The magnitude of Pq will 
be addressed later, but is not necessary for the following comparisons of different terms. 
The scalar potential can be cast in the form $o'^'r where 

47760 kpV kpT 

6 is the angle from the dipole axis. We can estimate $o/po ~ 3 • 10^^ J/C^m, while $r is 
again dimensionless. 

The Hamiltonian of the system is:— 

(-ffiV + eAr^_^^ (3) 

2m 



+ ^-^^ (4) 
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where e is the positive elementary charge. The A^-term can safely be neglected. Using the 
Lorentz gauge, V • A + c^<i> = 0, where c = (eo/io)^^^5 we rewrite the interaction Hamiltonian 

as 

//int = ArCa " V - ^.Cfe ; (5) 

Ca = 2/iBAo , Cb=(^l + e<l>o 

where /xs is the Bohr magneton. Since mc^ ^ /icu, C;, ~ e$o- 

For the wave function, we use the standard tight binding approximation,— writing the 
wave function as 

^k(r)=e^^--^EE^k.^nKr) (6) 

where \E'nz(r) = \E';(r — R„) is the l-th orbital corresponding to the atomic wave function 
centered at the ra-th atom located at Rn. The parameters 6ki can in principle be found for 
each point in /c-space. For states at the valence band maximum, there seems to be good 
agreement between theory and experiment. For the conduction band minimum, the fitting 
parameters are still optimized either for the position in fc-space or for the effective mass in 
different directions, depending on what is considered the most important. Based on Klimeck 
et a/.^'i^ we still assume that the minimum can be described by a combination of single 
electron p, s and s* states, where s* is an excited s-state. 

The transition rate for each k', k can then be found using Fermi's golden rule, 

Wk^k = — Kk'|i7i„t|k)|'5(^k' - - noo) (7) 

where k' denotes the final state, and k the initial one. In the following calculations, we will 
assume that the final state is near the conduction band minimum. There are six equivalent 
such minima, the effect of this will be addressed later. The absorbed energy by direct pair 
creation is given as Pd = J2k',k ^i^W^k',k- We calculate this using the parabolic approximation 
that is valid close to the band edges. Since we are interested in initial states close to the 
top of the valence band and final states close to a minimum in the conduction band one can 
assume that the interaction matrix element is weakly dependent on k and k', PVk'k ~ WkoO- 
Writing k" = k' — ko we get 

2mc ' 2m/ ^ ^ 



and using 



this gives 



Efe W^>^6{Ef -E,- fiu) (10) 



/ / -dEi^{Ef - Eg){-Ei)d{Ef -E,- huj) 

J En 



= W^„'^^§^(H.-E.f (11) 

where rric, rriy are the effective masses of the valence and conduction bands, respectively, and 
Eg is the gap energy. 

Note that the energy dependence of the absorption is the same as that of indirect absorp- 
tion, rather than that for direct absorption in direct band gap semiconductors, for which it 
is proportional to {ku; — EgY^'^. The reason for this is the spread of Fourier components in 
the dipole field, which take the role of the spread in phonon wave numbers in the case of 
indirect absorption. 

III. CALCULATIONS 

We term the power emitted as radiation Pr, the power lost to heating Ph, and the power 
going into direct electron- hole pair generation P^. 
After standard calculations we get the absorption: 
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2{mcmy)2uj{huj — Eg 
ttk'k = cos(|k' - k|r„) + ^ sin(|k' - k|r„) , 



1 (|k'-kp-fc^) + |k'-k| 



in 



Ik' - k| 



" I |k' -k\^-k^, + |k'-k|Va 1 • ^^^^ 

Here is the radius of the grain. {l'\0\l) = J 0^;(i^r denotes integration over the atomic 
orbitals for the operator O. It can be assumed that the elements of the sum where n' ^ n 
will only give small corrections. 
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Writing po ■ V = J2i=x,y,z^iPi'^ij and having an initial state that is a combination of p- 
states, only the matrix elements {s\eiPiV i\pi) do not vanish. They are expected to be of the 
order of 1/a where a is the lattice constant, 5.4 ■ 10~^° m for crystalline silicon. Regarding 
the contributions of the scalar field and Umklapp processes, only the elements with / = /' 
do not vanish. 

While each of the x, y, z give different contributions depending on the orientation of the 
dipole, it should be noted that there exist six equivalent minima in the conduction band. 
While the matrix element due to one minimum will be non-isotropic, the sum over all six 
minima is expected to be isotropic and equivalent to two minima with po || k. We define 
ko = |k' — k|. At the minima we have ^ 0.85 ■ 27r/a ~ lO^'^m"^. 

We see that all terms in Eq. [T^ show oscillations with period l/|k' — k| with increasing 
nanoparticle radius. As the nanoparticle diameter cannot be expected to be well defined on 
this length scale (atomic radius), we will simply take the average over one period. We believe 
this to be justified both from considering the limited coherence length of the electrons, and 
from the fact that any physical measurement would include a dispersion of particle sizes. 
While the limit of — > is mathematically well defined, it is not physically meaningful, as 
it describes a nanoparticle with less than one atom. 

Interestingly, the scalar potential provides much larger contribution than the vector po- 
tential. Keeping only the largest terms, we can make an order of magnitude estimate. 



Pd, Pr and Ph are all proportional to Pq. We define the damping coefficients 7^ = Pd/Po, 
■jr = Pr/Po, Ih = Ph/Po and 7 = 7d + 7r + Ih- We use these coefficients when comparing the 
importance of the different mechanisms. To find the total absorbed power, we also need to 
estimate the dipole moment, Pq, which is determined by the amplitude of the incident wave, 
Eq, and the polarizability of the nanoparticles, a, as po = a{uj)Eo. While the polarizability 
is in general dependent on the particle volume and shape, at the plasmon resonance it is 
determined by the damping only. This can be shown from equating the power absorbed 
from a plane wave by an oscillating dipole, EqPquj/2, with the total emitted power, 7Pq, 
giving 
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(13) 



00 



(14) 



a = 



27M' 
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Assuming that the incident hght is absorbed by a layer of nanoparticles with a 2D density 

n, we get 



Here W = eoElc/2 is the incident power per area. As long as 7^ is small compared to 7r + 7/i, 
it will not significantly change the polarizability, but if it becomes of the same order, the 
decrease in polarizability will be more important than the increase in absorption. This can 
easily be remedied by a higher density of nanoparticles, but if the density becomes very 
high, interaction between neighboring particles will change both the polarizability and the 
plasmon resonance frequency. 

IV. RESULTS 

We see that the direct absorption is proportional to {AEY = (^^^ — EgY) the excess 
energy after bridging the band gap, squared. This is the same energy dependence as found 
for indirect absorption if the single phonon processes including emission or absorption of a 
phonon are considered separately. 

The resonance of the nanoparticles has a certain width, so Ai? also has a spread. To get 
a feeling for the order of magnitude of Pd we define x = AE/hu and express the results 
through this. We are interested in frequencies where direct absorption of a plane wave 
would be impossible, so an x of close to one is irrelevant. The direct absorption shows no 
explicit dependence on temperature, as opposed to indirect absorption. This may indicate 
a method for differentiating between direct and indirect absorption in solar cells containing 
nanoparticles. The possibility of changes in the band structure of silicon with temperature 
should still be considered. 

The direct absorption, P^, should be compared with the energy lost to heating, Ph or by 
radiation, P^. From Eq. f|T2|) we have 



where x will usually be significantly less than one. With 50 nm particle radius we get 
kpTa ~ 0.3. The total integrated radiation from a dipole with the same dipole moment, po 




Id 



(15) 



oc 



{id + 7r + IhY ' 




(16) 
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Using that kpX = 27r we obtain 



^^"I27re " ^^^^ ' ^^^^ 



7. = ^ 10'' -4f-2 (18) 
127re sm^C 



giving 

^ ^ 3xM 1 + — ) . (19) 



7r V (kp^aY 



As shown in figure 1, we see that for sufficiently small particles and if x is not too small, 7^ 
can be of the same magnitude as 7,., or even larger. However, the gain for small particles 
requires that the particle is very close to the silicon. Designs where the nanoparticles are 
located outside the silicon may lose the benefit of the l/kpVa term. 

In the previously mentioned experiments, the nanoparticles were located close to a thin 
silicon wafer. For dipoles located at such an interface, the total radiation increases, and a 
large fraction of the radiation is directed into guided modes in the silicon wafer.-i^i^i^ For 
the light in the bound modes, it is assumed that the optical path length is sufficient to allow 
most of the radiation to be absorbed in indirect electron hole pair creations. At the same 
time the areas where the near field is strongest, has no silicon to absorb the energy. Under 
such circumstances, it should be assumed that 7^ dominates 7^. 

However, there is also a possibility that radiated energy can excite a plasmon on a neigh- 
boring particle, then again to be reemitted. This would reduce the positive contribution 
of the plasmons in architectures where the radiation along the wafer is exploited, as more 
energy would be lost to heat. If we instead place the nanoparticles inside the silicon and 
exploit the direct absorption, less energy would be reabsorbed by neighboring nanoparticles, 
and thereby less would be lost to heat. 

Using a simple resistivity argument, can be estimated as 

7^ ^ ^ 2 ■ 10^2^^—^—^ (20) 

for silver. For small particles, heating will take over as the dominant damping mechanism, 
as shown in figure 1. 
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V. CONCLUSIONS 



Our findings indicate that direct absorption due to surface plasmons on metal nanopar- 
ticles does occur, and may give important corrections to the total absorption for realistic 
parameters. 

The direct absorption has been found to be 

• proportional A£^^ = {hu — Eg)^; 

• independent of temperature; 

• inversely proportional to the k-space position of the conduction band gap minimum 
to the fourth power; 

• comparable in magnitude to radiated energy in some cases. 

The existence of the direct absorption mechanism is an argument for placing the nanopar- 
ticles inside the silicon, rather than in front of, or at the rear of the cell. This gives the 
additional requirement that the problem of recombination centers at the particle surface can 
be kept to a minimum. 

5 - 
4 - 

1044 

2 - 
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FIG. 1: The three damping mechanisms as function of kpr. For small kpr, 7/j dominates(dotted 
line), for large kpr, dominates (dashed line), while 7^ (solid line) may dominate in the middle 
region if x is sufficiently large (here for x = 0.15) 
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Ideally, the size of the nanoparticles should be so small that 7^ dominates 7^, but not so 
small that heating takes over as the dominant mechanism. We assume that diameters from 
about a tenth to half a wavelength could be suitable, depending on the conductivity of the 
nanoparticle. The plasmon resonance should be tuned using choice of material and particle 
shape (flattened for red-shift^), to match the band gap of silicon. 

The main questions that remain unanswered in our study concern the effects of interfaces 
and surface states for the direct absorption. We have not considered how an interface changes 
the near field, and the presence of surface states in the silicon may significantly change the 
problem in unpredictable ways. As the main contribution is from very near the dipole, 
both surface electron states and the alteration of the field due to an interface may be very 
important. There are also some unaddressed problems related to the averaging over particle 
radii and the finite coherence length of the electrons in the silicon. 

It is possible that higher order modes will give larger contributions to the absorption, 
these modes have been shown to be significant for nanoparticles of sizes where reradiation 
is larger than heating. 

The authors wish to acknowledge the financial support of the Norwegian Research Coun- 
cil. We wish to thank Alexander Ulyashin for introducing us to the idea of using nanoparticle 
surface plasmons in solar cells, and the SPREE at UNSW for giving us the necessary clues 
to get started. 



11 



^ p. Campbell, M. A. Green, Journal of Applied Physics 62, No. 1, Jul. 1, 1987, 243-249. 

2 P. Sheng, A. N. Bloch and R. S. Stepleman, Applied Physics Letters 43 (6), 579,(1983) 

3 K. R. Catchpole and S. PiUai, Journal of Applied Physics 100, 044504 (2006). 

^ H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings, Springer- Verlag 
(1988). 

^ C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles, John 

Wiley & Sons (1984). 

^ S. Pillai, K. R. Catchpole, T. Trupke, G. Zhang, J. Zhao and M. A. Green, Applied Physics 

Letters 88, 161102 (2006). 
^ S. Pillai, K. R. Catchpole, T. Trupke, M. A. Green, Journal of Applied Physics 101, 093105 

(2007). 

8 J. Mertz, J. Opt. Soc. Am. B 17, No. 11/November 2000 

9 H. Benisty, R. Stanley, M. Mayer, J. Opt. Soc. Am. A 15 No. 5/May 1998. 

° D. Derkacs, S. H. Lim, P. Matheu, W. Mar, E. T. Yu, Applied Physics Letters 89, 093103 
(2006). 

P. Lorrain and D. Carson, Electromagnetic Fields and Waves, Freeman 1970. 
^ L. D. Landau and E. M. Lifshitz Quantum Mechanics (Butterworth-Heinemann 2002). 
3 J. Singh Physics of Semiconductors and their Heterostructures McGraw-Hill (1993). 
^ G. Klimeck, R. C. Bowen, T. B. Boykin, C. Salazar-Lazaro, T. A. Cwik, A. Stoica, Superlattices 

and Microstructures, 27, 77 (2000). 
^ A. S. Martins, T. B. Boykin, G. Klimeck, B. Koiller, Phys. Rev. B 72 193204 (2005). 
^ J. Nelson The Physics of Solar Cells (Imperial College Press 2003). 
^ P. Wiirfel Physics of Solar Cells (Wiley- VCH 2005). 



12 



